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PREFACE 

This Memorandum is part of a series of continuing studies on 
orbital perturbations under NASA contract NASr-21 (02). It represents 
an extension to a previous RAND report, R- 399- NASA, "Perturbations of 
a Synchronous Satellite," by R. H. Frick and T. B. Garber. The 
results should be of interest to any agency concerned with long 
term satellite station- keeping. 




SUMMARY 


This Memorandum .analyzes the perturbations of circular equatorial 
orbits due to the ellipticity of the earth's equatorial section. The 
results indicate that with the exception of orbits with periods of 12, 
24, and 36 hours, the perturbations due to equatorial ellipticity are 
negligible and completely dominated by the perturbations resulting from 
initial condition errors. In particular, initial errors in orbital 
radius and orbital velocity cause steady state drifts in angular 
position of the satellite which are not influenced by the earth's 
equatorial ellipticity. 

The three special periods include the 24-hour synchronous orbit 
which, as shown previously, undergoes long period oscillations in 
longitude about the position of the minor axis of the earth's 
equatorial section. In addition, the orbits with 12 and 3 6 hour 
periods are resonant orbits in which the driving function resulting 
from equatorial ellipticity has a frequency equal to orbital 
frequency. This results in a divergent oscillatory perturbation at 
orbital frequency in both orbital radius and orbital angle. The rates 
of buildup of the amplitudes of these oscillations are 30.4 n mi/yr 
in radius and .243 deg/yr in angle for the 12-hour orbit. Hie 
corresponding figures for the 36 -hour orbit are 34.1 n mi/yr and 
.131 deg/yr.* These divergent oscillations are the dominant oscillatory 
effect for these two orbits, but they are still subject to the steady 
state drift in angle due to initial condition errors. 

*The exact values of these rates are dependent on the assumed 
ellipticity of the earth’s equatorial section. 
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SYMBOLS 



n 




o 

Ar 


Ar 

o 

Ar 

o 


6(Ar Q ) 

5(Ar Q ) 

t 

U 


6V 

o 

Y 


y 0 

Y 

9 

6 


gravitational acceleration at the earth's surface 
earth oblateness coefficient 

equatorial ellipticity coefficient 
modified, ellipticity coefficient 
ratio of orbital rate to earth rate 
mean earth radius 

radial distance from earth’s center to satellite 
circular orbital radius corrected for oblateness 
synchronous orbital radius 
circular orbital radius without oblateness 
perturbation in r 

amplitude of divergent oscillation in Ar 

initial value of Ar 

initial value of rate of change of Ar 

error in Ar 
o 

error in Ar 

o 

time 

earth' s gravitational potential 
error in initial orbital velocity 

angle between r and minor axis of earth's equatorial 
section 

initial value of y 
rate of change of y 

satellite central angle measured from inertial reference 
rate of change of 0 



X 


6 * earth rotation angle measured from inertial reference 

9 = earth's angular rate 

& 

= unperturbed orbital angular rate 

A0 = perturbation in 8 

= amplitude of divergent oscillation in LQ 

L8 - initial value of A 8 
o 

bQ = steady state value of A0 
ss 

6(A0 ) = error in A$ 
o o 
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I. INTRODUCTION 

It has been shown in Refs. 1 and 2 that a synchronous equatorial 
satellite under the influence of the earth’s equatorial ellipticity 
will undergo long period oscillations in longitude about the position 
of the minor axis of the earth’s equatorial section. 

Elis Memorandum considers the effect of this same ellipticity on 
circular equatorial orbits of any period. The analysis as presented 
in Section II is similar to the perturbation method used in Ref. 1 
with the generalization that the orbital angular rate no longer 
equals the earth’s rate. The solution of the perturbation equations 
shows the sensitivity of the resulting drift rate to initial condition 
errors. Three special cases of the general solution are also 
considered: the synchronous orbit and the two resonant orbits with 
angular rates of twice and two-thirds of earth rate. 
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II. ANALYSIS 


STATEMENT OF THE PROBLEM 

If a satellite is established in a circular equatorial orbit 
with an angular rate of 0 Q , what are the perturbations in angular 
rate and orbital radius resulting from the earth’s equatorial 
ellipticity? 


RkKkKEUCE SYSTEM 

In Fig. 1 the XY coordinate system is a geocentric inertial 
reference system in the earth's equatorial plane. The line AA' is 
the minor axis of the earth's equatorial section vhich rotates at 
earth rate, 0„. The instantaneous position of the satellite S is 

ill 

specified by its radial distance, r, and its central angle 0, 
measured from the X axis. The angle y is the instantaneous central 
angle between the satellite radius, r, and the minor axis AA .' . 


EQUATIONS OF MOTION 

In Ref. 1, the general equations of motion are developed. The 
potential function used is of the following form for the equatorial 
plane: 


U = 




1 + JL 


4 

2r 2 


+ 3J, 


( 2 ) ^ 


cos 2y 


( 1 ) 


The resulting equations of motion in the equatorial plane are of the 
form 


2 

dr *2 

— - re 

dt 


g o4 




9 J 2 <2) «o4 


2r 


cos 2y 


( 2 ) 
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Fig. I — Reference system 
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1 1 . ( r 2 b) = - 

r at Cr 0) 


^^4 


sin 2y 


(3) 


perturbation equations 

If the desired orbital angular rate is 0 Q , a set of perturbation 
variables can be defined as follows: 

r = r + Ar (*0 

c 

0 = 0 + A0 (5) 

o 


where r c is the circular orbital radius corresponding to b Q (r c is 
corrected for the effect of earth oblateness). 

Also the angle v can be expressed in the form 


V 


Y t + V 


where 



(6) 

(7) 


Substitution of Eqs. (4), (5), and (6) into Eqs. (2) and (3) gives 

* 

the desired perturbation equations as follows: 



4c (“> + !e (t 2 ' - & 2 ®o sin 2<Yt + Y o> 

c & o 


( 9 ) 


*For details see pp 6-8 in Ref. 1. 


where 
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k 2 = 


J ^ ** 

J 2 7 

c 


and the value of r c is obtained from the equation 

♦2 _ 

r c S o " — + ■■ 7 4 
r 2r 

c c 


The approximate solution of Eq. ( 11 ) is 


J 2 *E 

r = r l 1 + ■ ■ ■ g 

C ° V 2r 2 
o 


v ^ lere r o radius corresponding to 9 if there were no 


oblateness, as given by 


r 3 _ g°^E 


PERTURBATION SOLUTION 

Equations (8) and (9) can be solved simultaneously to give 

Ar 2n cos 2y n(5n-8) cos 2 y cos net 

— zz ^Ir 2 o E 


Ar _ „ “ w 

r ” ^2 n-1 

C L 


in- 2 ) ( 3 n- 2 ) 


2 n (n- 3 ) sin 2 y sin n 9 t n 2 (n- 3 ) cos 2 f(n-l) a t+v 1 
[n-2) (3n-2) Cn-l) (n-2) (3n-?) 


. Ar sin n e„t Ar 

+ (4 - 3 cos n e„t) — - + 5 £ 

E r * r 


c n 9 t , c 

ij 


A0 

+ 2(l-cos n o„t) — 7— 


( 1 ^) 
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e E 

2n (n-3) sin 2y Q sin n 6 ^ t 
+ (n-2) (3n-2) 

n (n 2 -2n+4) cos 2 [(n-1) + y q ] 

2(n-l) in-2 ) I3n-a) 


3n cos 2y n(5n-8) cos 2y Q cos n e fi t 
2(n-iT £ + (n-2) ( 3^-2 5 


• * 

Ar 2 sin n Q^t Ar 
„ * . \ O -C* __U 

6n (1-cos n 0„t) - ; ~ 

E r c e E 


A6 

+ (4 cos n flgt - 3) “T" 

e E 


(15) 


• • 

At * Ar and A0 arc the initial values of Ar> Ar^ and A0 
o 7 o o 

respectively* and n is given by 


n 



(16) 


An examination of Eqs. (l4) and (15) shows that they are 
indeterminate for n = 2/3, 1, and 2, corresponding to orbital periods 
of 3 6 , 24, a nd 12 hours. These cases are evaluated in the next 
section of the report. 

Ptor any other value of n than the three specified above, Eqs. 
(14) and (15) are valid solutions. By a suitable selection of 
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initial conditions, the coefficients of cos n Q„t, sin n 0_t and the 

E E 

constant term in Eq. (15) can be made zero. The required initial 
conditions are as follows: 


Ar Q 3k 2 n (n-3) cos 2y o 5(Ar Q ) 
r (n-l) (n-2) (3n-2) + r~ 


( 17 ) 


Ar 


6k 2 Q^ci (n-3) sin 2y Q 6(Ar ) 
(n-2) (3n-2) + r ° 


( 18 ) 


0 2 2 # 

A0 q 3k 2 n (n -2n+4) cos 2y 6(A0 ) 

T~ = (n-l) (n-2) (3n-2) 2 + 


E 


( 19 ) 


In each of these expressions the first term is the required initial 
condition while 5(Ar Q ), 5(Ar Q ) and 6(A0 q ) are the residual errors in 
achieving the desired value. Substitution of Eqs. (17), ( 18 ), and 
(19) into Eqs. (l4) and ( 15 ) reduces the solution to 

Ar _ ^2 1)2 005 2 [j 11 ” 1 ) + Vj 

~ r (n-l) (n-2) (3n-2) 


+ (4-3 cos n 0 t) 


. 8(Ar^) sin n Q^t S(Ar Q ) 


+ 


E r ■ r 

C n0 E 


5(A 0 ) 

+ 2 (l-cos n © E t) — — 


ne 


E 


( 20 ) 
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M 

0 


3k 2 n 2 (n 2 - 2n + 4) cos 2 [(n-1) 0 £ t + v Q ] 
(n-1) (n-2) (3n-2) 


E 


6n (l-cos n 0„t) 


8(ArJ 2 sin n 0 £ t 8(Ar Q ) 


E r 


0 


E 


5(A0 J 


+ (4 cos n 0 £ t - 3) — r 


( 21 ) 


An examination of Eqs . (20) and (21) shows that the earth's 

equatorial ellipticity causes a hounded oscillation with a frequency 

of 2(n-l) 0 in both Ar and A0. In addition the errors in the 

initial conditions introduce oscillatory terms at orbital frequency, 

n 0 as well as constant bias terms. In particular, the steady 
E 

state bias term in Eq. (21) would result in a steady drift of the 
satellite relative to the desired orbital rate, n 0 g. 

SPECIAL CASES 

It was indicated in the previous section that Eqs. (l4) and (15) 
become indeterminate for n = 2/3, 1 and 2. These cases are considered 
below. 

Synchronous Satellite (n = 1) 

If the limits of Eqs. (l4) and (15) are taken as n approaches 
1 the following expressions are obtained 
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— = 

r ^2 
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3 cos 2y - 3 cos 2y cos QJt 
o o E 


- 4 sin 2 y o sin 0g t + 4© E t sin 2y o 


• • 


. Ar sin e„t Ar 

+ (4-3 cos 0 t) — - + £ 2 

Ej r_ • r 


C 0 E 


A0 
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E 


( 22 ) 


A6 r 

— = ok 2 [- 3 COS 2V C + 3 cos 2 y ( 


o 003 V 


4 sin 2y o sin 0 £ t - 30 fi t sin 2 y q ] 


- 6 (1 


. Al% 

cos 0 t) 

E r 


2 sin 0 t Ar 

O 

• r 
0 E C 


A0. 


+ (4 cos e t - 3) -r- 2 


E 


(23) 


As before, the constant tern and the coefficients of cos 0 t and 
. E 

sin 0 £ t in Eq. (23) can be set equal to zero, so that the desired 
initial conditions are 


Ar o 

— = - 3k 2 cos 2 v q 


8 (Ar o ) 


r 

c 


(24) 
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4r . 6 < lr o> 

— - lac e E sin 2v 0 ♦ -y— 

c c 


( 25 ) 


6a 6(A0 J 

» 0 + — 


0 


E 


Or 


(26) 


These initial conditions reduce Eqs. (22) and (23) to 


~ = - 3k 2 cos 2 y q + 12k 2 0 E t sin 2y o 


+ (1+ - 3 cos 0 t) 


. 8(Ar ) sin 0„t 8(Ar ) 


E o 


6(A0 ) 

+ 2 (1 - cos 0 t) ; 

E *e 


( 27 ) 


A0 

e„ 


lSk 2 0 E t sin 2y o 


- 6 (l - cos 0„t) 


5(Ar Q ) 2 sin 0 £ t 8(Ar Q ) 

E " ~ r E ~ 


8(A0^) 


+ (4 cos 0 £ t - 3) — r 


E 


(28) 
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The residual bias term in Eq. (27) can be regarded as an additional 
correction to the steady state radius r,. Then it is seen that if 
the residual initial condition errors are reduced to zero, Eqs. (27) 
and (28) become 

r * 13t 2 V Sin (29) 

f = " l8k 2 V Sin ^o (30) 


These relations are equivalent to Eqs. (53) and (54) of Ref. 1. 


Twelve Hour Period Satellite (n = 2) 

If the limits of Eqs. (l4) and (15) are taken as n approaches 2, 
the following expressions are obtained 


~ = 3k |~4 cos 2y - 4 cos 2v cos 20_,t 
c - o o E 


+ \ sin 2y o sin 26 E t - e E t sin 2(0 fi t + y q )] 


. Ar sin 20 t Ar 
+ (4 - 3 cos 20 t) — - + 5 2. 

h T • r 

C 20 e r c 


A0 

+ (1 - cos 20 t) — - 


(31) 



and 

are 
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Ae 


= 12k [- 3 cos 2 y o + 3 cos 2 y o c ° s 


+ \ sin 2y 0 sin 20gt + 0 E t sin 2(0 E t + v Q )] 


Ar 2 sin 20 t Ar 

- 12(1 - cos 20 t) ; — 

E r = e E = 


A0 

+ (k cos 20 t - 3) ~ — 


(32) 


0 

As before the constant term and the coefficients of cos 20gt 

sin 20 t are set equal to zero. The required initial conditions 
E 

obtained, as follows: 


Ar 8(6r o ) 

-2 = - 3k 2 cos 2v 0 + — ^ 
c c 

Ar . 5 (Ar ) 

— " 3*2 S E sin + 

c c 


(33) 

( 34 ) 



o + 


5(A0 O ) 


(35) 


Substitution for the initial conditions in Eqs. (31) and (32) gives 
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Ax* • • 

— = - 3k 2 cx>s 2(e £ t + Y o ) - 3k 2 e E t sin 2(0^ + yj 


E T o' 


+ (4-3 cos 20 t) 


. . & ( Ar 0 ) sin 2 ® E t 6(Ar Q ) 


E r 


5(A0 1 ) 

+ (1 - cos 20„t) — 

tu 


20 


E 


0 


E 


(36) 


A0 


E 


12lC 2 e E t Sin 2 ^ Q E t + Y 0 ) 


12 (l - cos 20 t) 


. . 5 < Ar J 


2s in 20 t 5(Ar ) 
o 


E r 


8(A0 ) 

+ (4 cos 20 £ t - 3) — 

b E 


(37) 


From Eqs. (36) and (37) it is seen that both Ar and A0 contain 
oscillatory terms with amplitudes which increase with time. 

36 Hour Period Satellite (n = 2/3) 

If the limits of Eqs. (l4) and (15) are taken as n approaches 
2/3, the following expressions are obtained 
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- 4 cos 2v + 4 cos 2y cos 

1 ' O 


20 E t * 7 A 


+ ^ sin 2 y Q sin — ^ e R t sin 2\-y - Y, 


/ ' \ 20 E t * 

f „ 20 eM Ar o . 3 sin T Ar ° 

+ U - 3 cos — + : — 

\ 6 / c 20 e c 


+ 3 \ 1 - cos 


2e E t 1 A0 o 


= 4k 0 3 cos 2v - 3 cos 2y cos 


20 E t 7 A 


- ^ sin 2 y o sin — " 3 0 E t Sin 2 VT “ Y < 


- 4 l 1 - COS 


20 eM Ar o 2 sln 3 V Ar o 


3 / r 


f 20 E t \ A0 o 

+ U COS zr- - 3/ 

\ ' e E 


If as before the constant term and the coefficients of cos — =— 

2<Lt 

and. sin — j- in Eq. ( 39 ) are equated to zero, the following initial 
condition relations are obtained 

Ar 8(Ar ) 

— - + 3k_ cos 2v + — — M 

r 2 T o r 

c c 


ti . 8 (" 0 > 

— - - »£ e E sin + ~ r 
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Ae 

0 




E 


Qr 




Substitution of these initial conditions in Eqs. (38) and (39) gives 


Ar i E^ \ 

t = ^ ~ s s \— - v + ^ v eln 2 \rr - v 0 


2 Q *t 

/ 28„t\ S(4r ) 3 sin — 8(4r ) 
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(^3) 


A0 


28 * ( 

— 0 t sin 2* E 


3 2 E 
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28„t\ S(4r ) 


COS 


1] 

3 } r_ 


20„t 

2 sin — 5(Ar o ) 

• ~ 

0Tn C 


20„t \ 8(A0 ) 


+ Ucos -==- - 3 


0 
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( 44 ) 


As in the case of the 12-hour period satellite^ a divergent 
oscillation appears in both Ar and A0. 
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III. RESULTS AND DISCUSSION 


GENERAL SOLUTION 

If n is not equal to 2/3, 1 or 2, an examination of Eqs. (20) 
and (21) shows that the earth's equatorial ellipticity produces 
bounded oscillations in both Ar and A 6 at a frequency of 2(n-l) © E 
or twice the orbital rate relative to the earth. These oscillations 
result in displacements of the order of a few hundred feet relative 
to an unperturbed satellite with an orbital rate n0g. 

If in addition any residual initial condition errors exist, 
oscillatory terns will be introduced in both Ar and A 0 with a 
frequency of n0_, the orbital rate relative to inertial space. In 

hi 

addition, the initial condition errors will introduce constant bias 
errors in Ar and A Q which will cause a steady state drift relative to 
the unperturbed satellite. The magnitude of this drift is given by 

the relation 

5(Ar ) 

se ss = - 6ne E -j-2- - 38 <m o ) (*>5) 


where 



6(A0 o ) 


6V 


n 


2/3 


5V 


(W) 

m 


r = orbital radius for a synchronous 
s satellite (227^.^ n mi) 

5V = orbital velocity error 
o 
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Substitution of Eqs. (46) and (47) into (45) gives 

6e T 


A0 = n 5 ^ 3 &(Ar ) - 2 - n 2 / 3 5V 

ss r v o r o 

s s 


m 


Since the errors 8(Ar Q ) and 5V o are independent, the resulting 
contributions to the steady state drift can be treated separately 
so that 


60 


A0 


ss 


E n 5/3 


&(Ar Q ) 


= 34.66 n 


5/3 


6(Ar o ) deg/yr/n mi 




where 8(Ar Q ) is expressed in nautical miles. 

. 3 n 2 / 3 SV 

A6 

ss r 

s 

= - 39.19 n 2 / 3 8V 0 (50) 

where 8 V q is in ft/s. 

Figure 2 is a plot of Eqs. (49) and (50) as a function of n 
plotted from n = 0 to n = 17.067 which corresponds to a surface 
orbit with a period of 84.4 min. 

As an example, for n = 5, corresponding to an orbital period 
of 4 hrs and 48 min and a nominal orbital radius equal to 7779*9 
n mi, an error of .1 n mi in radius would result in a drift rate of 
50.7 deg/yr. 

Similarly an error of 1 ft/sec in orbital velocity would result 
in a drift rate of 114.6 deg/yr. 
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Fig. 2 — Steody stote drift rate 
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Hiis sensitivity of drift rate to initial condition errors is 
identical with that shown in Fig. 51 of Ref. 3* Hius, it is seen 
that the resulting drift rate is determined by the initial condition 
errors. Hie presence of equatorial ellipticity merely changes the 
desired values of the initial conditions as specified in Eqs. (17), 
(18), and (19). 

SPECIAL CASES 

In the special cases when n = 2/3, 1, or 2 it is seen that the 
solutions for Ar and A 0 include terms with amplitudes which increase 
with time. Hie physical significance of these time dependent 
solutions is discussed in the following paragraphs. 

Synchronous Satellite (n = 1 ) 

From Eqs. (29) and (30), it is seen that both Ar and A0 increase 
linearly with time. However, in Eq. (7) which defines y the 
additive term of A 0 has been neglected as being small compared to 
0 Q -0 E . In the case of a synchronous satellite, this is rc longer 
valid, since 0 Q is equal to 0^. If the A0 term is included in y, 
it effectively causes a slow change in y Q in Eqs. (29) and (30). 

Thus, these expressions for Ar and A0 are only valid for small values 
of A0 (of the order of 10°). However, in Ref. 1 a more general 
treatment of this problem is presented in which the limit on the 
size of A0 is removed. Hiis more general treatment shows that Eqs. 
(29) and (30) represent the initial part of a long period oscillation 
in satellite longitude about the position of the minor axis of the 
earth’s equatorial section. 




- 20 - 


Twelve- Hour Period Satellite (n = 2) 

This case and the one following for n = 2/3 are resonant orbits 
in which the frequency of the perturbing function, 2 y, is equal in 
magnitude to the orbital angular rate 0^« This condition is given 
by the relation 

2y = 2(n-l) e E = ± ne E = ± e c (51) 


which is satisfied for n = 2 or 2 / 3 . 

5br n = 2, Eqs. ( 36 ) and (37) show that both Ar and A0 contain 

terns of the fonn t sin 2 (© E t + y q ) which is the typical buildup of 

an undamped resonant system excited at its natural frequency. Thus, 

the amplitudes of Ar and A0 apparently grow without bound and the 

solution is only valid during the time that — and j— remain small. 

c 0 E 

It should be noted that the A0 contribution to y would detune 
the system slightly so that a beat frequency response would result 
between the frequencies 2y and 0 q . However, the initial buildup of 
this beat is represented quite well by Eqs. ( 36 ) and (37)* 

If the divergent tern in Eq. (37) is integrated, the following 
expression for the divergent term in A0 is obtained 

A0 = - 6 k 2 0 £ t cos 2(0 E t + v o ) (52) 


Thus for this orbit A0 has an oscillatory buildup at orbital 
frequency, 2 © E , and with an amplitude A 0 m which increases at the rate 



6k 2 S E 


= . 2^3 deg/yr 


(53) 
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wnere 



r 

c 



.6335 r 

5 


( 5 ^) 

(55) 


and 

= - 5.35 x 10 “ 6 (Ref. 4) 


(56) 


Similarly, it is seen from Eq. ( 36 ) that Ar also has an oscillatory 
buildup at orbital frequency with an amplitude Ar m which increases 
at a rate given by 



= 30.4 n mi/yr 


(57) 


If it is assumed that the initial errors in the radial and tangential 
velocities are both 1 ft/sec and that the initial error in orbital 
radius is 1 n mi, it can be shown that at the end of one year of 
operation, the divergent terns described above are the dominant 
oscillatory terms in both Ar and A© as expressed in Eqs. ( 36 ) and (37) 
However, the steady state drift rate in A © is still present, 
and its magnitude can be determined from Fig. 2 as 

A© = 110.0 deg/yr due to 5r 
ss J o 


A© 


ss 


= 62.2 deg/yr due to 5 V q 


and 
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Hiub a satellite with a 12-hour period would deviate from its unperturbed 
position with an oscillatory divergence superposed on a steady state 
drift rate in angular position. 


Thirty- Six Hour Period Satellite (n = 2/3 ) 

The behavior of the 36-hour period satellite is similar* to that 
described above. For the same initial position and velocity errors, 
it can be shown from Eqs. (43) and (44) that the dominant oscillatory 
terms in both fir and fi 0 are the divergent oscillations at orbital 
rate, 20^3, with amplitudes which increase at rates given by 


and 



lllk 2 K 


= 0.131 deg/yr 



* ^2 r c e E 


= 3^*1 n mi/yr 


(58) 


(59) 


As before, the steady state drift rate in A0 can be evaluated for 
the assumed initial conditions from Fig. 2 as 


and 


A0 =17.6 deg/yr due to 6r 

SS • Of* O 

A0 =29.9 deg/yr due to 8V 
ss o 


Thus, as in the case of the 12-hour period satellite, the deviation 
from the unperturbed position is an oscillatory divergence at orbital 
frequency superposed on a steady state drift rate in angular position 
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discussion 

Hie analysis presented in this Memorandum is similar to that 
developed by Blitzer in Ref. 5* in which he also investigates the 
perturbations of equatorial orbits, and in particular, the resonant 
12- hour and 36 -hour orbits. Hie rate of increase in the divergent 
amplitudes of Ar and Ad, determined in Eqs. (53)* (57)* (58)* and (59)* 
differ from those obtained by Blitzer, since the ellipticity which he 

-6 (2) 

assumed corresponds to a value of 1.67 x lcf for J^' ; based on Ref. 6 . 
If this value of is used in Eqs. (53)* (57)* (58), and (59)* the 

resulting numerical values are identical with those obtained by Blitzer. 

As a result of the analysis presented here and in Ref. 5* it is 
seen that for most equatorial orbits the perturbations due to equatorial 
ellipticity are negligible, the principal perturbations affecting 
station keeping are the steady state drifts in angular position due 
to initial condition errors. These steady state drifts are not 
influenced by the earth's equatorial ellipticity. 

The only cases in which equatorial ellipticity has any significant 
effect are the synchronous orbit and the two resonant orbits with 
12- and 36 -hour periods. The perturbation of the synchronous orbit 
is in the nature of a long period oscillation in longitude about the 
position of the minor axis of the earth's equatorial section and a 
small amplitude long period variation in orbital radius. 

Hie perturbations of the two resonant orbits in both angular 
position, Ad, and radial distance, Ar, are also oscillatory, but at 
orbital frequency and with an amplitude which increases slowly with 
time. Superposed on these divergent solutions are the usual 
perturbations resulting from initial condition errors. 
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It should be noted that the orbital rates have been specified 
relative to inertial space. Thus, the orbital rate as seen from the 
rotating earth is given by (n-l) 9„. In the case of the 12-hour 

ill 

orbit (n = 2) , the relative orbital rate is + 0 and the satellite 

ill 

moves toward the east and is over the same equatorial position 
every 24 hours. Similarly for the 36-hour orbit, n = 2/3, the 
relative orbital rate is - b^J 3 and the satellite apparently moves 
to the west and is over the same equatorial position every 72 hours. 
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IV. CONCLUSIONS 


As a result of this analysis, the following specific conclusions 
can be stated: 

o Hie only significant effect of the earth's equatorial 

ellipticity on equatorial orbits is for orbital periods of 
12, 24, and 36 hours. 

o The 2k - hour or synchronous orbit has been discussed in detail 
in Ref. 1. 

o The 12- and 36 -hour orbits are resonant orbits which display 
a divergent oscillatory perturbation at orbital frequency, 
o The amplitude of the divergent oscillatory perturbation in 
angular position grows at the rate of .243 deg/yr for the 
12-hour orbit and .131 deg/yr for the 36 -hour orbit, 
o The amplitude of the divergent oscillatory perturbation in 
orbital radius grows at the rate of 30.4 n mi/yr for the 
12 -hour orbit and 34.1 n mi/yr for the 36 -hour orbit, 
o For any orbital period other than 12 , 24, or 36 hours, the 
perturbations due to equatorial ellipticity are negligible, 
o From the point of view of long term station-keeping, the 
drift rates due to initial errors in orbital radius and 
orbital velocity are the dominant effect and may result in 
large displacements in angle from the unperturbed position. 
These steady state drifts are not influenced by the earth's 
equatorial ellipticity. 
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